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With the improvement of the experimental techniques, many new phenomena in which the photon 
degrees of freedom are involved have been discovered. Typical examples are the exciton-polariton 
quasi-particles, excitons strongly coupled to photons. A correct description of these systems re¬ 
quires a full account of the photon’s dynamics, including the vector gauge degrees of freedom. In 
order to include this contribution in ab initio many-body perturbation theories, here we present a 
generalization of the Hedin equations, originally derived for a many electron system with Coulomb 
interaction. These equations are now derived to take into account the strong coupling physics of 
the electrons and the transverse photon’s degrees of freedom. 

PACS numbers: 


I. INTRODUCTION 


In the past decade, progress in designing semiconduc¬ 
tor microcavities has attracted much interest from differ¬ 
ent fields of physics, ranging from condensed matter to 
quantum optics and relativistic particle physicJ^^. Typ¬ 
ically, a microcavity is a system where a two dimensional 
(2D) quantum well (QW) is embedded in Bragg mir¬ 
ror barriers. The confinement of the photons enhances 
the external electromagnetic (EM) interaction with the 
QW so that retardation effects due to the EM vector po¬ 
tential cannot be considered negligible, even though the 
system is in a non relativistic regime. When strongly 
coupled to the photons, elementary boson excitations 
can give rise to new quasiparticles: polaritons. Typi¬ 
cal examples are the exciton-polaritons, photons dressed 
by electron-hole excitations bearing a very light effec¬ 
tive mass (iTie ~ 10“^). Other polaritons stemming from 
the ultra-strong coupling regim^ were recently the sub¬ 
ject of intensive investigations due to the possibility of 
tailoring their properties by tuning the size of the QW, 
and its potential use for laser polariton device^^. The 
ultra-strong coupling regime has been observed in doped 
semiconductor QWs, where plasmons originated in the 
interconduction subband re sult in stable and robust in¬ 
tersubband (ISB) polaritond^l^lfiSI A similar effect takes 
place in the presence of an applied magnetic field, where 
the photon mode is coupled to 2D electron gas magneto- 
plasmon excitationJSI. Since polaritons are themselves 
bosons, Bose-Einstein Condensation (BEC) phenomena 
can be observed and invest igatecP. 

From a theoretical point of view, microcavity quantum 
electrodynamics (Q ED) syst ems are usually described 
by analytical model^^^Ultlll, whereas first-principle ap¬ 
proaches have been proposed only recently, extending 
the time-dependent density-functional theory (TD-DFT) 
to the relativistic regime^^. Any concrete applications 
of TD-DFT methods still require an appropriate ap¬ 
proximation scheme for the electron-photon exchange- 
correlation (xc) functionalJi^. 


Recently, an xc potential for such systems has been 
proposed within the time dependent optimized effec¬ 
tive potential theory (TD-OEP) frameworl^i^, enabling 
the possibility of first-principles calculations. In con¬ 
trast, ab initio many body perturbation theory (MBPT) 
approache^i^, which together with TD-DFT, have been 
very successful in describing the excited states of con¬ 
densed matter systems, have not been used so far to study 
the microcavity problem. Ah initio MBPT calculations 
are based on the self consistent solution of the Hedin 
equations, originally developed for studying electron in¬ 
teractions meditated by the instantaneous Coulomb in¬ 
teraction. These equations were then generalized to take 
into account spin-orbit and spin-spin interaction^^ or to 
be implemented within the dynamical mean field theory 
framework (DMFT)!^. 

In this paper, we present the generalization of the Hedin 
equations to non relativistic systems of charged particles, 
where the quantum degrees of freedom of the transverse 
EM gauge field is taken into account. The new Hedin 
equations, which are solved self-consistently, are the ba¬ 
sis of possible ab initio MBPT calculations, allowing to 
obtain the electron and photon quasi particle self energy. 


In section we introduce the interactions in the sys¬ 
tem following the gauge principl^23flll^ as done in QED. 
From the gauged Lagrangian, and in order to make con¬ 
nection with the condensed matter literature, we then 
obtain an effective Hamiltonian solely defined in terms 
of the fermionic degrees of freedom. In the second part, 
we use the effective fermionic description to formulate the 
problem in terms of the Schwinger fu ncti onal approach 
discussed in Ref.^^. Finally in section HI we apply our 
formalism to the historical case of the exciton-polariton 
in microcavity QED. Conclusions are drawn in section 
|Vj The formal derivation of the Hedin equation is pre¬ 
sented in the appendix. 
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II. THEORETICAL DERIVATION 


The Lagrangian density for a system of non relativistic 
electrons is our starting point; by demanding invariance 
with respect to C/(l) gauge transformations this reads 


^ _ jeao) + -leAi)'^ 




( 1 ) 


Here V'(l) is a fermionic (Grassman) held, ao(a:) is the 
internal scalar potential (due to electron-electron in¬ 
teractions) and Ai is the external (laser) EM gauge 
held. We dehne the total EM gauge held as A^{x) = 
{c ao{x), Ai{x)) and the bare g auge held propagator as 
9 *^ ( 1 ^ 2126 ] j^g^Ye also introduced 
the hat metric tensor rj^i, and the index 1 = (a;i,ti,s) 
stands for space, time and spin variables. In the follow¬ 
ing we will use Einstein’s convention for summation over 
dummy indices and /r = (0, i), i = x, y, z (the temporal 
and spatial components, respectively) and units of fi = 1. 
In our reference frame, x and y are the transverse direc¬ 
tion and z is the longitudinal direction. Nevertheless, as 
we focus on non-relativistic systems, the covariant and 
contravariant indices are indistinguishable. Choosing the 
Coulomb gauge V • A = 0, the propagator of the gauge 
helds is diagonal, and in momentum space it reads^S 
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( 2 ) 


where one can recognize the Coulomb propagator in the 
top left corner. Equation Q can now be written in the 
equivalent form 


— !> V’ + poo-f jM, (3) 

- {^) 

where p(x) = jo{x)/c = e V'^(x)'!/'(x) is the charge density 
and ji{x) = 1/(2toz) [tpHx) (Viz/>(x)) - (Viz/;'l'(x)) ipix)] 
is the (paramagnetic) charge current. Note that contrary 
to the relativistic theory, a term proportional to (dia¬ 
magnetic term) now appears. The effect of this term is 
to renormalize the spatial component of the gauge field 
propagator. It is then convenient to define the renormal¬ 
ized gauge field propagator as 

(-^(0))/2!^ “ (^(Q) )/jy H (4) 
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FIG. 1: Diagrams contributing to the photon propa¬ 
gator: the first terms in the diagrammatic expansion of the 
renormalized photon propagator, Eq. 


Hmb(1) =H°(1)+?(1)C^(1) 

-II d2fil)iD^o))nAh2)r{2), 


(5) 


where "H® is the non-interacting Hamiltonian density 


-H°(l) = z^1'(l) + ^(’'i)| V’(l) 


( 6 ) 


In Eqs. ([^ and (|^ we have introduced the crystalline 
nuclear potential V and the source field needed to gen¬ 
erate current expectation values. Finally, note that we 
have moved to the operator formalism and we will desig¬ 
nate operator valued quantities with a hat. In the effec¬ 
tive description, the absence of the gauge field is traded 
for the non local current-current interaction term. If the 
spatial components of the current are neglected, the well 
known result for the instantaneous Coulomb interaction 
in terms of a non lo cal density-density interaction term 
is recovered | 22 | 24 | 26 | 

The one and two point electron Green’s functions are 
defined in the same way as the time ordered expectation 
value. 


Gi(l,2) = (iV 

T 

z^(l)i/)1'(2) 

v) 

G2(l,2;T,2') = (iV 

T 

z/i(l)V)(2)V;’l'(2')z/’'*'(l') 


where \N) is the many body ground state. Given the 
above definitions, the charge and current densities can 
be expressed as 


p(l) = -zeGi(l,l+) 


( 8 ) 


a(i) 



Gi(l,l') 


e;(i,i 0 Gi(i,i'), 


where the propagator should be understood in terms of 
its diagrammatic expansion (see Fig. 0). Integrating 
out the gauge fields, we arrive at the effective Hamil¬ 
tonian density containing only the fermionic degrees of 
freedom. 


where 1+ stands for -|-£ (e —)• 0+). In order to be con¬ 
sistent with this notation ^o(l) — C(5(l, 1'). From the 

Hamiltonian density of Eq. ([^, the equation of motion 
(EOM) for the electron Green’s functions is obtained in 
analogy to those in Ref.^^ 
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Gi( 1 , 2 )- J d3 (^Gi(l,3)Co(3)(5(3,2) + |iGi(l,3)Cz(3)5(3,2)^ 

+ i I d3(^p(o))oo(l,3)G2(l,3+;2,3++) + (Z?(o))„ (l,3)2l,(3+,3+)G2(l,3+;2,3++)^ =^(1,2). 


This EOM is comprised of the zero-component part de¬ 
scribing the usual longitudinal scalar part (T>(o))oo (the 
Coulomb potential) and the new transverse components 
due to the EM gauge field. The formal derivation of 
Hedin equations is shown in the Appendix. Here we re¬ 
port the final results: The Dyson equation for the one- 
particle electron Green’s function has the usual form: 

Gi(1,2) = g(°>(1,2)+ f d(34)Gl°^(l,3)E(3,4) Gi(4,2), 

( 10 ) 

where E(l, 2 ) is the self energy operator and G^^^ is the 
Green’s function for the non-interacting Hamiltonian. 
The Dyson equation for the self-energy E is 

E(l,2) = 

-i J d(3456)r(o)^(6,5;l)D^,(6,3)Gi(5,4)r"(4,2;3). 

( 11 ) 

Note that in these equations, we have introduced the 
dressed Photon propagator (as before Dgo is the 

screened potential usually written as W) and T"" is the 
irreducible vertex operator . Therefore, the ”GW” self 
energy is now retrieved by replacing the vertex T'^ with 
the bare one, defined as 

r(0)M(i, 2; 3) = 4(3, 3')5(1, 3')5(2,3). (12) 

The Dyson equation for the irreducible vertex operator 
is 

r'^(l,2;3) =r(°)'^(l,2;3)-t 

/ ^|^^Gi(4,6)Gi(7,5)r'^(6,7;3). 

(13) 

The Dyson equation of the screened photon Green’s func¬ 
tion is 

l^/ii'(l,2) =(Z)(o));i,y(l, 2)4- 

/ d(34) (D(o));.^(l, 3) n^e(3,4) Dfl,(4, 2), 

(14) 

whereas the Dyson equation for the polarization function 
is 

n^.(i,2) = 

f d(4567) r(°)'^(7,6; 1) Gi(6,4) Gi(5,7) r'^(4,5; 2). 

(15) 


The set of Eqs. pl)|-( 15) (summarized in Table. 1) con¬ 
stitutes the non relativistic Hedin’s equations in which 
the full photonic degrees of freedom are taken into ac¬ 
count. Some further remarks are needed here: 


1. In this framework, the Bethe-Salpeter equa¬ 


tions ( [l3| and (151 coincide with Eqs. (C.ll) and 
(C.29) already derived in Ref.^^. 


2. The /ormaZ structure of these non-relativistic QED 
Hedin’s equations is analogous to two cases al¬ 
ready studied in the literature: The first is the 
relativistic cas^^, derived by replacing the vec¬ 
tor current with the relativistic charge current 
ju = '|f(a;) 7 ^ 4>(a:) (where now 4>(a;) are spinors 
and 7 ^ are the Dirac 7 -matrices). As we have al¬ 
ready explained, in this case there is no diamag¬ 
netic contribution however. Einally, in the spin 
dependent cas^^, one needs to replace the para¬ 
magnetic current with the magnetization operator 
m* = 'I'l(a:) (T® ^(a;) (with cr* being the Pauli ma¬ 
trices). 


3. The ’’diamagnetic” term in the photon Green’s 
function Dq contributes new diagrams (see Eig. [^, 
similar to the quantum back reaction in scalar 

qeeESI. 



FIG. 2: semiconductor microcavity Scheme. The QW 

is embedded in two Bragg mirrors The wave vector along z 
direction is quantized whereas the is free. 
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TABLE I: The Hedin equation in the presence of the external transverse EM field 


Summary table of the Hedin equations 
Gi = Gf ^ S Gi 


pM = p(0)M + ^GiGiT'" 

= r(°)'^GiGir" 


III. THE EXCITON-POLARITON 


In a QW embedded by two Bragg mirrors (microcav¬ 
ity), the wave vector in the z direction perpendicular 
to the cavity plane is quantized, whereas the in-plane 
motion is free (see Figj^. This brings us to an impor¬ 
tant physical phenomenon: the coupling between a pla¬ 
nar cavity photon mode and a boson such as an exciton, 
which leads to a quasi-particle called a polariton. This 
makes an ideal case in which the generalized Hedin equa¬ 
tion can be applied. In order to describe the exciton- 
polariton, we, first start from the free transverse photon 
Green’s function for the photons generated by the exter¬ 
nal source and trapped in the microcavity. In reciprocal 
space, Eq. ([^ reads 




4:Tr[Sij - {kjkj/k'^)] 
uj^ — + irj 


(16) 


where 77 0 is a regularization parameter. From the 

Hedin equation (141 and by using the expression for 
ZI(o) ill Et- (161) the dressed transverse photon Green’s 
function ia^^I 


- (kjkj/P)] 

— c^k"^ + ir]’ 


(17) 


where we have introduced the macroscopic, transverse 
dielectric function 





47re^ 


Hfj (fe, uj). 


(18) 


Here the plasma frequency Up = is ob¬ 

tained from the diamagnetic correction in the bare pho¬ 
ton Green’s function by keeping only the first order term, 
see Fig[l] The transverse polarization 11 can be inter¬ 
preted as the self energy of the photon in analogy to the 
standard case. The dispersion relation of the exciton- 
polaritons are obtained from the poles of D in Eq. (|17||^ 

det |a;^e"^(k, w) — = 0. (19) 

Considering the z direction as the axis perpendicular to 
the QW, the photon motion along this direction is quan¬ 
tized as kz = ttM/Iz, with M being a positive integer 


and Iz being the thickness of the cavitji^. Therefore, in 
Eq. (19) k'^ = {ttM/IzY + ^jj- initio calculations of 
dispersion relations for excitons have recently been per¬ 
formed by solving the Bethe Salpeter equation, show ing 
that nowadays they are computationally feasibld2213iJ_ 


IV. CONCLUSIONS 

In conclusion, we have generalized the Hedin equations 
in order to take into account the dynamics of the EM 
transverse field for systems in the non-relativistic regime. 
In order to do it, we have obtained an effective Hamil¬ 
tonian in which the degrees of freedom of the internal 
EM gauge field are integrated out. This leads to a new 
term coming from the diamagnetic current in the photon 
Green’s function. By using the Schwinger approach, the 
generalized Hedin equations are derived. The extension 
proposed here is general in nature. As a possible applica¬ 
tion, we have focused our attention on the self-consistent 
solution of these equations in the context of resonant cav¬ 
ities, where the photon-matter coupling is enhanced by 
photon confinement. In this regard, we have considered 
one of the most important phenomena taking place in 
QED microcavities: the exciton-polariton, for which a 
many body Green function has been obtained. Neverthe¬ 
less, because our formalism is completely general, it can 
also be used for ab initio studies of polariton-polaritorP^ 
and photon-photon effect^. The ab initio study of the 
ultra-strong coupling regime of ISB polaritons will be the 
focus of future investigations. 
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Appendix A: Derivation of the Hedin equations in 
microcavities 

In Eq. Q the interacting terms (both temporal and 
spatial terms) are a direct sum in the space-time space. 


Therefore, they can be addressed separately. Since the 
temporal component leads to the known Hedin equations, 
in this Appendix we focus only on the new terms stem¬ 
ming from the external vector field Ai. Recalling that 
the vector source field is (i, we have 


^;(3, 3)G2(1, 3; 2,3+) = Gi(l, 2)C.(3,3)Gi(3,3+) - 

oQ(3) 

Here the mass operator M is defined as follows: 

M(l,3)Gi(3,2)= yd3(G(o))„(3,l)24(3+,3+)G2(l,3++;2,3+)= J Gi(l, 2)j,(3)- 

= I d3S(l,3)Gi(3,2) 

where we have introduced the self energy S and is the Hartree term: 

V/^{1) = J d2(G(o)),,(l,2)j,(2) 

With the total field: , the inverse dielectric tensor is 

-Hi 9 l = ^Mi)-A An -A An 914 . f 

= <5(1,2)+ J d34(G(o)),fe(l,3)nfc,(3,4)e,-i(4,2) 


(Al) 


<5Gi(1,2) 

<50(3) 


(A2) 


(A3) 


where we have also defined the polarization nfe;(l, 2) = <5jfc(l)/<5<i);(2). The Hedin photon Green’s function, Eq. (141, 
is obtained by using Eq. (A3) and considering that 


A,(1<2) = I <i3e-0l,3)(G(o))fe,(3,2) 

A, (1,2) = (^(o))^,(l,2) + I d345(G(o)).fc(l,5)nfc,(5,4)e-^i(4,3)(G(o)),„,(3,2) 

The irreducible vertex is defined as 


3*. (3) 


and with the chain functional derivative, we have 

^Gr 0 i, 2 ) 


A,(3) =y^4r(l,2;4)e-(4,3) 


(A4) 

(AS) 

(A6) 

(A7) 


The electronic self-energy E Dyson equation © is obtained from Eq. ( |A2[ ) (considering that the dressed photon 
Green’s function is defined in Eq. (A4): 


E(l,2) = -* I d35(D(o)).,(3,l)|^^^^Gr05,2) 


2m <50(3) 

then by using the identity (<5Gi/<50 = —Gi (<5G("^/<5C) Gi, we find 

<5Gr04,6), 


E(l,2) = -* J d3456(D(o)).,(3,l)|^Gi(l,4)^^^y^Gi(6,5)Gr05,2) 


(AS) 


(A9) 
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and we finally arrive at 


E(l,2) . -/«(%))«(3.1)gG.(I,4)ia^. 


(AlO) 


The Hedin equation ( fTTI) is recovered from Eq. (AlO) by using the bare vertex definition equation (12), Eq. (|A7[ ), and 
the definition in Eq. ( |A6[ ). In this framework, the Hedin vertex equation is the same as that shown in Ref.l^^l 

r (1,2; 3) = 6(3, 3')<5(1, 3')<5(2,3) + / d4567Gi(4,6)Gi(7,5)rX6,7; 3) = 


<5Gi(4,5) 


= r(o)Xl,2;3) + I d4567^^^Gi(4,6)Gi(7,5)r(6,7;3) 


(All) 


Similarly, the Dyson equation for the polarization is derived, considering the definition = 

6Ul)/6<^,{2): 


n,, (1,2) =6(1,1') 


^Gi(l,T) 

< 5 $, ( 2 ) 


d346(l,l')Gi(l,3) 


S^,{2) 


Gi(4,l') 


(A12) 
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